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Abstract. The Krauss-model is a stochastic model for traffic flow which is continuous in space. For periodic 
boundary conditions it is well understood and known to display a non-unique flow-density relation (funda- 
mental diagram) for certain densities. In many applications, however, the behaviour under open boundary 
conditions plays a crucial role. In contrast to all models investigated so far, the high flow states of the 
Krauss-model are not metastable, but also stable. Nevertheless we find that the current in open systems 
obeys an extremal principle introduced for the case of simpler discrete models. The phase diagram of the 
open system will be completely determined by the fundamental diagram of the periodic system through 
this principle. In order to allow the investigation of the whole state space of the Krauss-model, appropri- 
ate strategies for the injection of cars into the system are needed. Two methods solving this problem are 
discussed and the boundary-induced phase transitions for both methods are studied. We also suggest a 
supplementary rule for the extremal principle to account for cases where not all the possible bulk states 
are generated by the chosen boundary conditions. 

PACS. 02.50.Ey Stochastic processes; 45.70.Vn Granular models of complex systems, traffic flow; 05. 40. -a 
Fluctuation phenomena, random processes, noise, and Brownian motion 



1 Krauss-model 

The number of vehicles on highways and in cities is in- 
creasing each year causing vehicular traffic to suffer more 
and more from jams. The phenomena related to traffic 
jams have attracted the attention of physicists and engi- 
neers since almost half a century, trying to develop models 
describing the features of the real traffic. Generally there 
are two different approaches: microscopic and macroscopic 
||l|,|[H]. Whereas in microscopic models different vehicles 
and their dynamics can be distinguished, in macroscopic 
models only densities are considered, similar to hydrody- 
namcis. 

However, the approach of a physicist is usually quite 
different from that of a traffic engineer. One of the current 
interests of statistical physicists are the so called " nonequi- 
librium systems" . In microscopic vehicular traffic theories, 
vehicular traffic is treated as a system of interacting parti- 
cles driven far from equilibrium and offers the possibility 
to study various fundamental aspects of the dynamics of 
truly nonequilibrium systems. 

Empirical observations show that the average velocity 
decreases with increasing vehicle density. So the average 
current (or flow), which is the product of average velocity 
and density, is a function of the density. The functional re- 
lation between current and density is usually called funda- 



mental diagram. Its generic form can be understood easily. 
For small densities all vehicles can move with their desired 
velocity Wmax a-iid the current increases monotonously. For 
large densities the vehicles interact with each other and 
the average velocity is much smaller than in the free flow 
regime. This causes a decrease of the current, with a max- 
imum at an intermediate value. 

The traffic model introduced in [^|5[^ , called Krauss- 
model in the following, is based on an approach by Gipps 
[0 considering the braking distance of individual cars. 
Starting from the assumption of safe driving an update 
scheme can be formulated in the manner of the well-known 
Nagel-Schreckenberg (NaSch) model [|[||]. In the Krauss- 
model — unlike the NaSch model — the state variables, 
i.e. space and velocity, are chosen to be continuous. To 
make the model safe, i.e. free of collisions, a safe velocity 
Vsafc for each car is introduced, which is calculated in ev- 
ery timestep taking into account that there is a maximum 
acceleration and deceleration rate for each car. The vehi- 
cles will be updated in parallel corresponding to discrete 
time dynamics. 

The model has been designed to reproduce the empirial 



findings in traffic jams [hO 
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1. There is a density regime with non- unique flow-density 
relation. 
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Fig. 1. Classes of qualitatively different behaviour in the 
Krauss-model. In class /, the acceleration a and deceleration 
b are realistic and all properties f — 4 of jams are present. In 
class II decelerations are large and the properties 1,3, and 4 
of jams are not reproduced. In class /// accelerations are large 
and no jams exist (from |5|). 



— the stability of the outflow interface of jams. 

These stability criteria were introduced to obtain a clearer 
characterization of traffic flow models with the focus on 
their stochastic properties. The Krauss-model of class / 
exhibits stable high-flow states and a stable jam interface 
[ p^ . Stability of high-flow states means that the intrinsic 
stochasticity of the model is not sufflcient to cause a tran- 
sition into the jammed regime. Therefore, the dynamics 
of the Krauss-models differs from the VDR- model |p^ , ^ , 
for which an unstable interface was found. For the latter, 
the high-flow states are truly metastable, i.e. for increas- 
ing system length the probability for a transition into the 
jammed state becomes equal to one. In order to emphasize 
the difference between the nature of the high-flow states 
in the two models the term bistable is used in the context 
of the Krauss-model. 



1.1 Dynamical equations 



2. Traffic jams can develop and exist under "pure" con- 
ditions, i.e. in the absence of any obstacles. 

3. The flux out of a jam is not maximal. 

4. The outflow from jams is stable. 

5. The outflow from jams and the velocity of the down- 
stream front do not depend on the inflow conditions. 

These properties are displayed by the model for a certain 
range of parameters. The model equations proposed in 
even show a much richer behaviour depending on the re- 
lation between the ac- and deceleration capabilities . 
Three different domains can be distinguished (see Fig. m: 

— Class /: 

— Accelerations and decelerations are realistic and 
bounded. 

— All properties of jams are modeled correctly. 

— The jamming transition is a first order phase tran- 
sition. 

— The interactions are effectively long ranged. 

- Class //: 

— Decelerations are unbounded, leading to effectively 
short ranged interactions. 

— Properties 1, 3, and 4 of jams are not reproduced. 

— The jamming transition is no phase transition, but 
a crossover. 

- Class ///: 

— Accelerations are unbounded. 

— No structure formation at all. 

Throughout this article only models of class / will be in- 
vestigated, i.e. stable jams can occur and the flow-density 
relation is not unique in a certain density regime (Fig. 
The stability of the jams is directly related to the fact 
that the outflow from a jam is smaller than the maximal 
possible flow |l3|| . 

Recently an alternative classification of stochastic traf- 
fic models with respect to two properties has been sug- 
gested 

— the stability of the high-flow states. 



To derive the underlying dynamical equations, two types 
of motion of vehicles are considered. The first type is free 
motion, the second the motion of a vehicle while interac- 
tion with another vehicle takes place. Corresponding to 
this, two main assumptions can be made. The free motion 
is bounded by some maximum velocity Wmax^ 

V < Wmax ■ (1) 

It is assumed that the system remains free of collisions 
and that a driver always chooses a velocity that does not 
exceed the maximum safe velocity Wsafo which guarantees 
the absence of collisions: 

V < Wsafo ■ (2) 

fsafc is determined from the condition that the braking 
distance d{v) needed to stop when moving with velocity v 
satisfies 

d{vf) + VfT < d{vi) + g . (3) 

The quantity on the left side is the braking distance of 
the following car (velocity Vf) including a finite reaction 
time T. This distance has to be smaller than the braking 
distance of the leading car (moving with velocity vi) plus 
the gap g between the vehicles. Furthermore the model 
takes into account that positive and negative accelerations 
are bounded: 

-b < — < a , with a, 6 > 0. (4) 

It is natural to implement these restrictions using a con- 
tinuous space variable, but time is discrete with timesteps 
At. From the above restrictions the dynamics of the model 
can be derived. 

It will be assumed that, apart from random fiuctua- 
tions, every vehicle moves at the highest velocity compat- 
ible with the restrictions stated above. In this way the 
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model can be formulated immediately, giving 



2b 

Vdcs{t) = min{i;max, v{t) - 
v{t + At) = max{0, Wdcs(i) - v} , 
x{t + At) = x{t) + vAt . 



+ T 

- aAt,Vss,{c{t)} , 



(5) 



Here the gap g — xi — Xf — Zcar is the spatial headway 
between the leading car at xi and the following car at a;/ , 
where ^car denotes the length of a car. vi and w/ are the 
velocities of leading and following cars, respectively. The 
safe velocity f safe has to be determined in accordance with 
condition (0) . v^cs is the desired velocity representing the 
wish to drive as fast as possible through the acceleration 
V + a At, but also respecting the conditions (|l|) and (^. 
The random perturbation 77 > has been introduced to 
allow for deviations from optimal driving, where rj = 
and ^ is a random number uniformly distributed in the 
interval [0, 1]. 

In the following we set At = t = \. The other generic 
parameters used in the simulations are 



a = 0.1, b = 0.6, Wmax = 5, e = 1.0, ^car = 1-0 



(6) 



The unit of the space coordinates is the length ^car of 
one car. Another parameter is the length L of the sys- 
tem which has been choosen to be equal to 2001 (if not 
stated otherwise). 



1.2 Characteristics of the model 

For the parameters chosen in (^) the Krauss-model be- 
longs to class I and exhibits a bistable region with a stable 
high flow branch which implies a non-unique flow-density 
relation. Fig. ^ shows a fundamental diagram for a system 
corresponding to class /. 

The existence of a bistable regime is related to the oc- 
curance of phase separation in the system. The distribu- 
tion of the gaps and velocities of a system in the jammed 
state has two peaks i.e. there are two groups of cars 
in the system. It separates into a macroscopic jam and a 
free-flow region. According to the initial conditions there 
exists another system state in which cars drive with veloc- 
ities close to Wmax and the distribution of gaps possesses 
only one peak. These states belong to the high-flow branch 
in the ambiguous part of the fundamental diagram. 



Fig. 2. Fundamental diagram of the Krauss-model with pe- 
riodic boundary conditions and a — 0.1, b = 0.6, Wmax = 5, 
e — 1.0, Zcar = 1.0. The density is a mean value of densities 
measured in an interval of length located in the middle of a 
system of length L = 2001. The flow q is measured locally in 
the middle of the system. For densities 0.1 < p < 0.14 the flow 
is bistable such that a stable high flow branch can be observed. 



parameters, namely the injection rate a and the extrac- 
tion rate f3. So the car density in the bulk will be a re- 
sult of these rates and the underlying model dynamics. In 
general a nontrivial density profile will develop, i.e. the 
average density in the system will depend on the position. 

The influence of a and /3 on the car density implies 
that quantities like bulk density^, current (flow) q, and 
the density profiles show a different behaviour than in pe- 
riodic systems, which were studied extensively for cellular 
automata, optimal velocity models etc. (see e.g. [|l|,|[|^ 
and references therein). 

On the other hand, for OBC most investigations deal 
with simple one-component systems, especially the asym- 
metric simple exclusion process (ASEP) |l7, 18, 19|,p0|,pi|. 
In the NaSch model with Wmax > 1 was studied with 
OBC and the results compared to empirical data. Special 
boundary conditions for this case were also studied in 
0. For driven lattice systems which exhibit a metastable 
or bistable regime for periodic boundary conditions not 
much is known about possible phase diagrams in the case 
of open boundaries (see, however, |16,2j]), nor have sys- 
tems been studied with stable high flow branches. 



2 Open boundary conditions 

One of the most significant differences between systems 
with open and periodic boundary conditions is the car 
density p, which in a periodic system is a conserved quan- 
tity. Here the density and the initial conditions (in the 
bistable regime) determine the stationary state completely, 
which allows to study the density-dependence of the macro- 
scopic parameters. In systems with open boundary condi- 
tions (OBC) one has to deal with two different tuning 



2.1 ASEP with open boundary conditions 

The asymmetric simple exclusion process (ASEP) is the 
simplest prototype-model of interacting systems driven far 
from equilibrium. It is a generic model for studying driven 
systems and boundary-induced phase transitions p5|,[2q. 



^ Measured by averaging over an interval of length ^ in the 
middle of a system of length L. 
^ See, however, the discussion of these results in |l^. 
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The ASEP is a discrete particle hopping model. A par- 
ticle can move forward one cell with probability p if the 
lattice site immediately in front of it is empty. If the first 
cell, corresponding to the left boundary^, is empty a par- 
ticle will be injected there with probability a. If the last 
cell is occupied the particle will be removed with proba- 
bility /3. By varying the tuning parameters a and /?, and 
therefore the densities at the boundaries, one obtains a 
surprisingly rich phase diagram. 

One distinguishes three different phases according to 
the functional dependence of the current and the corre- 
sponding stationary bulk density on the system parame- 
ters. In the low-density phase the current is independent 
of p. Here the current is limited by the input rate a which 
then dominates the behaviour of the system. In the high- 
density phase the behaviour is dominated by the output 
rate /? and the current is independent of a. In the maxi- 
mum current phase the limiting factor for the current is 
the bulk rate p and the current becomes independent of 
both a and /?. 

In ||2^ a nice physical picture has been developed which 
explains the structure of the phase diagram not only qual- 
itatively, but also quantitatively. By considering the col- 
lective velocity Vc — q'{p) which is the velocity of the 
center of mass of a local perturbation in a homogeneous, 
stationary background of density p and the shock velocity 
of a 'domain wall' between two stationary re- 



_ g2-gl 



P2-P1 



gions of densities pi and p2, one can understand the phase 
diagram of systems with unique flow-density relation from 
the fundamental diagram of the periodic system ||2^,p9{. 
The idea behind is that these two velocities determine if 
and how a perturbation will spread through the system. 
For a detailed discussion see |29[ . 

A general valid "rule" is found for systems with unique 
flow-density relations, i.e. the current always obeys an ex- 



tremal current principle ||28,|29 



g = maxpg[p+_p-] for p > p^ 

q = minpg[p-,p+] q{p) for p- < p+ 



(7) 



p~ and p"*" are effective densities at the left and right 
boundary, respectively. The principle (|^) states that the 
phase diagram of the open system is completely deter- 
mined by the fundamental diagram q[p) of the periodic 
counterpart. Moreover, it implies that two models with 
different microscopic dynamics, but the same fundamen- 
tal diagram, will have the same phase diagram for open 
boundaries. In this sense the phase diagram is indepen- 
dent of the microscopic dynamics. 



2.2 Krauss-model with open boundary conditions 

The aim of this paper is the study of the Krauss-model 
with open boundary conditions, especially obtaining its 
phase diagram, and furthermore to investigate its connec- 
tion to the theory of boundary-induced phase transitions 



(see section 2.1). As will be shown, the choice of appropri- 
ate inj ection/ extraction strategies at the boundaries of the 
system plays a crucial role. Here it should be kept in mind 
that the principle @) is formulated in terms of effective 
boundary densities which result from these strategies. 

The rules specified in the following sections have to be 
such that the full range of possible bulk states (compare 
Fig. H) may be reached (at least theoretically). Especially, 
we are interested in states of high flow. Therefore, one 
might think of a strategy to inject cars with a initial ve- 
locity of Wmax and injection rate a. However, this will not 
lead to a crash-free motion by itself. 

Since the model is known to be crash-free from the 
closed system, this is somehow surprising. The reason can 
be found from the fact that under open boundary condi- 
tions all kind of initial situations can occur due to the 
stochastic feeding of cars. For the safety of the model 
quantity ^(i) = g{t) — vi{t) plays a cruicial role. Its evo- 
lution for the deterministic Krauss-model (e = 0) is given 

by § 



i{t + At)>m 1 



1 



n + i 



(8) 



We assume that the particles move from left to right. 



with T}, = [v -\- vi)/2b. Equation (^) implies that if once 
^ > (and as a result g > vi > 0) this will hold for all 
future timesteps. Safety is therefore guaranteed ii £,{t = 
0) > 0. The latter condition is not fulflUed automatically 
if cars are fed with a rate a. In simulations we found that 
a car which collides with its predecessor always had ^{t — 
0) < 0. Note that the opposite is not true, i.e. a car that 
started with negative ^ does not have to be involved in 
a crash ^. Due to the stochastic step in the update rules 
(||), ^ can be pushed from negative to positive values but 
not the other way around. Just choosing smaller initial 
velocities reduces the probability that £,{t = 0) < leads 
to a crash (even to negligible values) but the states of high 
flow will not be reached ^. 

We follow two different strategies of injecting cars to 
overcome the mentioned problems. One strategy is based 
on the following idea: If there is enough space at the be- 
ginning of the system, i.e. at least one carlength ^car, cars 
are injected according to the injection rate a. As a conse- 
quence one has to define an initial velocity Winit = Vf(t = 
0) that is as high as possible (since in the high flow states 
the average velocity v « Wmax), but keeps the system free 
of crashes. The problem here is, when using the formula 
for fsafc from (|^), that a velocity Vf is already needed for 
its calculation. Moreover one has to deal with the cases in 
which £,{t — 0) < 0. Using the safety condition (||) we will 
derive a rule to determine ■^init in section |[ 

The other strategy investigated goes the opposite way. 
The high-current states of the Krauss-model are charac- 
terized by velocities close to fmax and more or less iden- 
tical gaps g ^ Icar- This property can be used to define a 

* From (^) follows that a crash is most probable for slow 
moving leaders and small initial gap ginit. 

^ That is another difference to the NaSch-like models in 
which the choice of v < g always leads to collision-free mo- 
tion. 
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rule which mimics the structure of the high current states, 
i.e. one defines a minimal gap ginit > har that has to be 
respected at the left boundary and injects all cars with 
t^init = timax- It should be noted that in this case the injec- 
tion rate does not equal a anymore, but becomes a mono- 
toneously increasing function of a. Details and simulation 
results are given in section H. 



3 An inflow-oriented injection rule 

In this section an injection method is introduced which is 
similar to that for the ASEP. Cars are injected into the 
system with inflow rate a whenever there is at least one 
carlength space in front of the system {gmit > ^car), using 
any safe initial velocity ?;init (depending on the system 
configuration). However, in order to reach large currents, 
cars have to be injected with the maximum safe velocity 
possible. 

The rule, as stated up to now, leads the condition 
^{t = 0) = ginit < which can cause accidents as seen 



in section 2.2. Since we want to investigate a naive gen- 



eralisation of the injection strategy of the ASEP and to 
compare the results to it, it is necessary to think about 
a rule which is mainly oriented on a. Moreover, in real- 
world applications one has to understand the behaviour 
of such a rule, since cars usually are inserted according to 
a given inflow instead of particular strategies. 



3.1 Boundary rules {rule 1) 

In order to complete the rule we have to find the maximum 
safe velocity Winit possible. We can not just use (||) since 
Vf is not known. Moreover, the velocity has to be such 
that the dynamics of the system allows the transition from 
6nit < to ^ > 0. As a solution we do not look only at 
the first car in the system, but also at its predecessor. 

The open boundary conditions for the inflow-oriented 
rule are defined in the following way: 

Step 1: Injection 

If there is at least one car length free space at the begin- 
ning of the system, with probability a we inject a car with 
velocity Winit: 



■^init 



mm < v-n 



2bg 



(9) 



This velocity is a function of Wmaxj the velocity vi of the 
leading car, the deceleration rate b and an upper bound 
bi for the actual deceleration of the leading car. The lat- 
ter is calculated using the velocity of the car in front of 
the leading car and therefore, gives a bound for the worst 
case, i.e., the maximum deceleration of the leader in the 
next timestep. 
Step 2: First update 

Performing the first update of a car injected at the cur- 
rent timestep, we define an own rule. Given Vinit of step 1 
we follow the update rules of the Krauss-model in case 



of the leading car moving with vi > Wcrit, where Wcrit is a 
constant velocity depending on Wmax- If ^ ^critj ^safe is 
set equal to the gap g instead of using (||). This defines a 
cutoff for which ^init < still leads to safe driving while 
keeping ^init close to 1 for high values of a. 

Step 3: Extraction 

With probability 1 — /? a block is added at the end of the 
road which causes the car at the end of the system to 
slow down. Otherwise, with probability (3, the cars simply 
move out of the system. 

Step 4- Update 

Update with the Krauss-model update rules (see section^. 

In the following we use parameters as given in (|^) and 
fcrit = 1-6, .ginit = ^car = 1-0. The valuc for Wcrit has been 
determind by means of simulation. Note that ginit is the 
space that has to be free at least at the left border of the 
system. 



3.2 Fundamental diagram 

In contrast to periodic boundary conditions, the funda- 
mental diagram q(p) is not easy to find for the full range 
of bulk densities. While for a closed system the density 
p is given and conserved, in open system it is a quantity 
that results from the parameters a and (3. Their influence 
on p or g is non-linear. Another difficulty is that global 
density and current should be measured in the stationary 
state which is reached after quite long simulation times 
for certain values of (a,/3). A detailed discussion can be 
found in Because of the complex relation between 
(a, (3) pairs and p or one can not find a value for the q 
for each p and vice versa. It should be mentioned that p 
and q are rather sensitive to changes in a or f] [^o| . 

Our strategy in finding the fundamental diagram is as 
follows: For each pair of parameters (a, /?), the simulation 
has to be run until the stationary state is reached. Then, 
the flow q and bulk density p are measured in the middle 
of the system. For the latter this has been done in an 
interval of length . These two values fix a point in the 
p-q plane. To make the data more reliable and to reduce 
the influence of noise, this has been repeated for several 
times with different random seeds. 

In Fig. 1^ the fundamental diagram of a system with 
periodic boundary conditions and the same system with 
open boundary conditions are compared. The high-flow 
branch in the bistable density regime does not exist in 
the open case. There are also other states missing in the 
fundamental diagram of the system with open boundary 
conditions, i.e., certain densities can not be generated by 
the boundary rules defined in section 3.1. 



3.3 Phase diagram 

The phase diagram as a function of a and /3 can be ob- 
tained by studying the density profiles of systems in the 
stationary state. It is easy to distinguish between low- and 
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Fig. 3. Comparison between the fundamental diagram of the 
open system (■) and the one from the closed system (o). The 
points are obtained from density profiles in case of the open 
system. 
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high-density phases. To find the maximum current phase 
one should compare the (p, q) pair of the system with the 
fundamental diagram of the same system with periodic 
boundary conditions. 

Using the density profiles for systems with a, /3 G [0, 1] 
with increments of 0.02, we have drawn a phase diagram in 
Fig. ^. Different phases, phase boundaries and the typical 
density profiles for each phase are shown. The maximum 
current phase with q about 0.5 has been observed for sys- 
tem with open right boundary (/? — 1) and injection rate 
between a = 0.48 and a = 0.57 (broken line). Note that 
all stable high-flow states in the bistable region correspond 
to some point in the maximal current phase. 

This phase diagram looks different from phase dia- 
grams of similar models since one finds two different low 
density regimes. Therefore, by varying a with (3 kept con- 
stant, one observes a reentrance transition for large values 
of /3. It is natural to assume that this is related to the spe- 
cial choice of input and output strategies. To verify this, 
later on (Sec. ^ a different injection strategy will be stud- 
ied. The reason why this phase diagram looks unusual can 
be understood studying the relations betw een a, P and 
dynamic parameters of the system (see Sec. 3.4). 

A similar phenomenon has been observed in | |3l[ | in a 
simple one-component lattice gas with next-nearest neigh- 
bour interaction. Here a reentrance transition to a second 
high-density phase was found. Although the origin of this 
transition is not entirely clear, the authors of pi[ argued 
that it is related to the complicated connection between 
boundary rates and the effective boundary densities. 



3.4 Dependence of dynamic quantities on a, (3 

The dynamics of a system with open boundary conditions 
depends on the parameters a and (3. When /? = 1, i.e. 
the outflow is unrestricted, an increasing flow can be ex- 
pected with increasing injection rate a. In Fig. ^ the flow 



Fig. 4. Phase diagram for the Krauss-model with open bound- 
ary conditions and the inflow-oriented injection strategy. The 
inserts show typical density profiles for each phase. The max- 
imum current phase with flow q ~ 0.5 is only observed for 
systems with open right boundary (/? = 1) and injection rate 
between a = 0.48 and a = 0.57 (broken line) and shows an 
oscillating density profile. 



of the stationary state is plotted versus a. However, this 
is true only for a < 0.58. For larger values, a sharp de- 
cline in the flow can be seen which can be related to the 
injection strategy. Note that the second injection strategy 
introduced in section ^ does not show that decline in the 
a — g-relation. 

For systems with (3^1 the right boundary introduces 
an external disturbance which increases with decreasing 
(3 (cf. step 3). With (3^1, the cars at the end of the 
system are forced to slow down. Hence, the density at the 
right boundary increases with decreasing (3 and jams are 
formed. These jams grow backwards into the system. If 
the jams can not dissolve due to high inflow rates a, the 
cars have a lower average velocity and the system's density 
is high. 

Because of the bistability, the system has quite a dif- 
ferent behaviour for /3 = 1 in comparison to any other 
(3. Note that already small (external) perturbations might 
force the breakdown of the stable high-flow states in the 
open system. Once a jam has established in the system 
the high-flow states will not be reached again due to the 
reduced outflow from jams, (p, gout) = (0.11,0.51), which 
is smaller than the maximal possible flow. As an example 
for the behaviour of the system under a weak disturbance, 
the a — g-relationship for j3 = 0.95 is shown in Fig. ^. 

On the other hand, one can study the relationship be- 
tween q and (3, using a constant value of a (not shown). 
It is obvious that for very small values of /3, the current 
should vanish. In a system with closed right boundary 
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Fig. 5. Of — g-diagram for /3 G [0.70, 1]. Tlie maximal reached 
current decreases for smaller /3 since jams are builded due to 
the disturbance at at the right boundary of the system. The 
strong drop of the current for a > 0.58 is a result of the cutoff 
Wcrit in step 2. 



Fig. 7. Simulation o f th e left boundary rules given in sec- 



tion 3.1 




Fig. 6. Average velocity as a function of a and (3. The contour 

line ( ), given hy v > Umax/2, separates the free flow phase 

(average v close to Umax) from the jammed phase {v much 
smaller then Wmax/2). 



(Rule 1) and 4.1 (Rule 2). In each timestep only the 
leftmost three cars are simulated while the rest of the system is 
cutted off. The dots represent the fundamental diagram taken 
from the closed system. After a relaxation of 100000 timesteps 
the points (p, q) are measured over another 100000 steps for 
different values of a. The arrow indicates the direction of in- 
creasing a. 



or velocity, i.e., the rest of the system is cutted off. The 
resulting relation between the density and flow at the left 
boundary is shown in Fig. ^ Following the line starting 
from a = 0.1 the left border stays on the free flow branch 
of the fundamental diagram up to a = 0.6. For bigger 
values the system switches to just one state that belongs 
to the jammed branch. The drop in the flow is therefore 
not alone a result of jams moving backwards to the left 
boundary, but an artiflcial effect of the rule itself. If the 
predecessing vehicles are moving slow {v < Ucrit) f\ the 
safe velocity of the inserted cars becomes g due to step 2 
of rule 1. And, g will be very small in case of high injection 
rates. These effect will play a role for the interpretation of 
the phase diagram in the context of the extremal principle. 



3.5 Extremal principle 



(/3 = 0), all cars are forced to stay in the system which 
means a vanishing average velocity v = Q and density 
p = 1. For each a, one expects the highest value of cur- 
rent for /3 = 1. 

In Fig. ^the average velocity as a function of a and (3 is 
shown. The contour in the a-/3-plane shows the line where 
V = Wmax/2. In 1^ this line was choosen to distinguish 
between jam and free flow, i.e. states with v < Wmax/2 are 
in the jammed phase. 

Before we will examine the connection to the extremal 
principle the effects of rule 1 are investigated in more de- 
tail. In order to measure the impact of that rule on flow 
and density we let the system run for different values of a. 
In each timestep only three cars are left in the system by 
taking out the rest without taking care on their position 



Since we could not observe the full range of states (cf. 
Fig. ^ due to the cutoff in rule J at a certain inflow we 



examined the extremal principle (see Sec. 2.1) only for a 
subset of system states. Since the extremal principle (|^) 
is formulated in terms of boundary densities these have to 
be determinded from the stationary density profiles. For 
small injection rates it is not easy at all to find the left 
boundary density p~ , due to the oscillations in density 
profiles. This is no problem for models that have Umax = 
^car {— particle size). 

For the parts of the phase diagram that look similar 
to other models, i.e., all except for the second low density 
regime a > 0.6, /3 > 0.6 (cf. Fig. it should just be 

^ which results e.g. from jam formation around the left 
boundary 
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Fig. 8. Density profile for the system with (a, B) — (0.8, 0.8) 
which lies in the second low density phase of Fig. W. The inserts 
show a close up at the boundaries. The densities at the left 
and right boundary are p~ — 0.65 and = 0.f7 respectively, 
the bulk density is pb ~ 0.04. A stationary flow q = 0.21 is 
measured. 



stated that the principle was fulfilled for all these pairs of 
(a,/3). ^ 

A different behaviour is found for the second low den- 



sity regime. The results of section 3.4 suggest that this 
is due to the fact that for rule 1 we do not have a mono- 
tonic increase of inflow with a. Consider the point (a, (3) = 
(0.8, 0.8) which lies in the second low density phase. Fig. || 
shows the corresponding stationary density profile. Using 
the formulation (|^) one finds an apparent violation of the 
extremal principle, as will be demonstrated. From the den- 
sity profiles we find = 0.65 and = 0.17 and since 

> p+ the maximum of q{p~) and q{p^) will be chosen 
according to (Q), i.e., ^(p™'^) = 0.48. Instead we measure 
a stationary fiow oi q — 0.21. 

Before we present an alternative interpretation of the 
results in the second low density phase we recall some 
findings in the context of the ASEP. All the results have 
been obtained for models with Wmax = 1 = At. There- 
fore, in these models the flow obtruded to the system by 
the boundary rules does play no role. A stopped particle 
can accelerate to Wmax = 1 in one timestep. Then, the 
formulation as given in (Q) is sufficient to determine the 
system state. For models with Umax > 1 (cf. the left 
boundary rule has to be defined in a specific way to ob- 
tain results compatible with the extremal principle. The 
rule used in Jl6t always allows the injection of cars with 
'^max- Hence large flows can be reached and the problem 
that acceleration to the maximum velocity takes several 
timesteps does not occur. 

In our case the left boundary density is large due to 
the high injection rate {p~ oc a) but the flow is restricted 
io q = 0.21 due to the cutoff in step 2 of rule 1 (cf. Fig. ^. 
Therefore, the system has to choose a state that matches 
with that flow. Since the exit allows a higher flow no stable 
growing jam can develope in the system and the state on 
the free flow branch is choosen. Indeed the profile shows 



a sharp decrease in the density leading to a bulk density 
Pb = 0.04 for which q{p^,) = 0.21 (cf. Fig. |). This in- 
dicates that the extremal principle ([7|) formulated only in 
terms of boundary densities p", p^ is not sufficient. More- 
over, one has to check if there is a restriction due to the 
infiow and outfiow rules (denoted by g"''"'"). (0) then only 
applies for q < q~l^ . Otherwise the system state is chosen 

by 



qh 



i{q-(a),(7+(/3)} 



(10) 



with a density p^ satisfying q(^ph) ~ %■ Predicting the sta- 
tionary system state according to (^ together with ( p!c| ) 
one finds good agreement for all pairs of (a,/?). The re- 
sults given for another injection strategy as formulated 
in section ^ will confirm this interpretation. We believe 
that similar extensions of the extremal principle will be 
necessary for other multicomponent models. 



4 High-velocity-oriented injection rule 

As seen in the last section it is not easy to define left 
boundary conditions that allow the system to reach states 
belonging to the stable high-flow branch and - at the same 
time - ensure safe motion of the injected car under any 
circumstances. In the following a rule is formulated that 
tries to inject cars with the maximal velocity Umax of the 
model. 



4.1 Boundary rules {rule 2) 

In order to generate the high-flow states we have a closer 
look at the bistable regime of the corresponding peri- 
odic system. Here high-flow states are characterized by 
velocities close to Vmax and approximately identical gaps 
g 3> /car- Therefore, instead of driving with initial speed 
Winit to achieve safety, whenever there is an empty space of 
one carlength /car in front of the system, we try to inject 
cars (with probability a) with an initial velocity Umax- To 
guarantee a system free of crashes, we have to introduce 
a minimum safety gap ffi nit t o the preceding car (cf. also 
the discussion in section pT^ ). Therefore not all injection 
trials will be successful (see below). 

Step 1 and step 2 of rule 1 are replaced by above strat- 
egy while the deflnition of the right boundary conditions 
is not changed. Using this strategy, it might happen that 
no car will be injected for several timesteps because of 
the lack of free space (ginit) at the beginning of the sys- 
tem. This implies that the actual injection rate is smaller 
than a and also the existence of a maximal density — 
which can be reached at the left boundary. 

For the choice of the initial gap ijinit three criteria are 
formulated: 

(i) 5init should leave the system crash-free. 

{ii) High currents (comparable to those found in the bistable 

regime of the periodic system) should be reached. 
{Hi) (7i„it should be as small as possible. 
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Fig. 9. The relation between the current q and a for several 
values ginit G [1, 3.5], /9 = 1 and the parameters equal to (0). 



The last criterion increases the maximum reachable 
density at the system's entry so that a larger range of 
boundary densities can be investigated. 

To find an appropriate value for the initial gap, we 
made several simulations with different values for 5init, 
and measured the current q for (3 — Ito test criterion {ii). 
For the values of ^init which satisfy {ii) we then checked 
in simulations whether the system is free of crashes. This 
could be done best using /3 « 0. After filtering out unsuit- 
able values of g-mit we chose the minimum of the remaining 
values. 

In Fig. ^ the relation between q and a for several val- 
ues of ginit and /3 = 1 is shown. The current increases as a 
increases and either reaches a constant value or decreases 
and then reaches a constant value, depending on ginit- For 
5 > (?init > 2 the current reaches a constant value greater 
than 0.6, which is in the bistable region of a system with 
periodic boundary conditions. For ijinit > 5 q decreases 
drastically. This is obvious, since it will be impossible to 
reach a density in the regime of maximum current in that 
case. The maximum current in terms of the periodic sys- 
tem will be reached for a system with (3=1 and density 
of 0.15. 

In simulations for systems with [3=1 and /? <C 1 
we have observed crashes for g^^it < 2. Therefore, the 
only values of gi^a meeting the criteria (i), (ii) are val- 
ues 5 > (jinit > 2. Because of the criterion (Hi) we decided 
to take ginit = 2 in order to reach a maximum possible left 
boundary density. 

In the following we present results of simulations using 
this injection rule, parameters as given in ^ and ginit — 
2, Winit = Wmax- As long as the methods and interpretations 
correspond to section a the presentation will be kept brief. 




Fig. 10. Fundamental diagram of the system with open bound- 
aries and boundary conditions defined by rule 2. The maxi- 
mum current and the high-flow states in the bistable region 
are reached for systems with /? = 1. The fundamental diagram 
is similar to the fundamental diagram of the system with pe- 
riodic boundary conditions (see section hi). 



using rule 2, see Fig. |T^. With the present injection strat- 
egy, increasing a does not cause the decrease of the current 
in the way we saw in section ^. The high-flow states are 
only reached for systems with (3=1 (cf. next subsections). 



4.3 Dependence of dynamic quantities on a, (3 

In this section we give a brief overview of the dependence 
of the current q, average velocity v and the density p, 
measured in the middle of the system, on (a,/3). 

In Fig. |ll| the relation between q, a and /3 is shown. 
The current increases strongly for /3 values close to 1. We 
have also indicated the boundary between the regions with 
q > 0.5 and q < 0.5. 

Fig. |l2| shows the dependence of the density on (a, /3). 
The density used here is measured in the middle of a sys- 
tem in the stationary state, which is the density corre- 
sponding to the plateau in the density profile for most 
pairs (ct,/?). One cleary sees a sharp transition from low 
to high values which will be interpreted as the phase tran- 
sition line between high and low density phases. 

Finally in Fig. |l^ the dependence of the average ve- 
locity on (a, /?) is presented. Using again the criterion 
V = t'max/2 to distinguish between free flow and lam one 
finds the same results for the phases as in Fig. |l2. 



4.2 Fundamental diagram 

The fundamental diagram obtained for this system is very 
similar to the system with periodic boundary conditions. 
The high-flow states in the bistable regime can be reached 



4.4 Phase diagram 

In Fig. |lj the phase diagram and different phases are 
shown. The diagram has been derived from the density 
profiles. The full Q;~/?-plane was scanned in steps of size 
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Fig. 11. The relation between the q, a and (3 (ginit = 2). There 
is no decrease of the current with increasing a anymore. The 
broken line correponds to g = 0.5. 




Fig. 12. The relation between the density p (measured in the 
middle of the system in stationary state), a and /3 (ginit = 2). 
The broken line, corresponding to a density p — 0.15, can be 
interpreted as a phase transition line between high and low 
density phases. 

0.1. The maximum current phase has been observed for 
systems with open right boundary (/? — 1). It is reached 
for a > 0.5. As one can see from Fig. R at this value 
cars are fed in the system with the flow out of jam, i.e., 
(/in(a = 0.5) = 0.51 = gout 

This phase diagram is similar to that of the much sim- 
pler driven lattice gas discussed in |Q. The results com- 
pare as well to the phase diagram found for the VDR- 
model ||T^. Therefore, we will only briefly discuss the ma- 
jor differences found. 

For the driven lattice gas of the maximum current 
phase is only found for very short system. The reason is 
the true metastable nature of the high-flow states, i.e. in- 
trinsic fluctuations are able to destroy these states even 
without an external disturbance. Since the probability for 

^ Using a = 0.5 and /3 G [0, 1] one obtains the full high 
density branch of the fundamental diagram. 




Fig. 13. The relation between the average velocity, a and /3 
(ffinit = 2). The broken line corresponds to « = Wmax/2 and can 
be interpreted as transition line between free flow and jammed 
phase. 




a 

Fig. 14. Phase diagram of the Krauss-model with open bound- 
ary conditions and maximum-velocity-oriented injection rule 
{rule 2). The maximum current phase has been observed for 
systems with open right boundary (/3 = 1) and a > 0.5. The 
inserts show (from left to right) typical density profiles for low, 
high and maximum current phase. 

such fluctuations grow with the system length the max- 
imum current phase will vanish above a typical system 
size. In our case, we do not find such a disappearance, 
even for very large systems since the high-flow states of 
the Krauss-model are stable (as found in |0). 

The fact that we flnd the maximum current phase only 
for f3 — 1, while in the case of the VDR- model it exists 
for a slightly bigger range of /3, can be related to the fol- 
lowing. The interaction in the VDR-model is very short- 

* We checked this for systems up to L = 50000. 
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ranged due to the unbounded deceleration capability in 
that model (the interaction horizon is g ~ Wmax)- Instead, 
for the Krauss-modcl of the investigated class /, cars do 
already interact for g w u^^x- 



4.5 Extremal principle 

We examined the validity of the extremal principle for this 
model for every a and f3 G [0, 1] with 0.1 steps. After mea- 
suring the boundary densities and we have used 
to determine the current g(P''^'^) predicted by the extremal 
principle. This value has then been compared to the cur- 
rent q measured in the open system. For all parameters 
the results are in excellent agreement with g(P™^. This is 
in accordance with our suggestion (see section |3.5D that as 
long as the boundary rules do not restrict the flow in an 
artificial way, the system will choose its bulk state accord- 
ing to the extremal principle as formulated in (^ . We also 
checked our formulation with the addition given in (p^. 
For all pairs of (a, (3) the correct bulk state was predicted. 



5 Conclusions 

In this article we investigated the properties of a traf- 
fic flow model introduced by Krauss et al. under 
open boundary conditions. It is a discrete map in time 
while — in contrast to cellular automata approaches — 
space is continuous. For the parameter range discussed, 
the Krauss-model shows a non-unique relation between 
flow and density. This property of the fundamental di- 
agram is responsible for the stability of jams found in 
empirical observations. It is important to note that the 
corresponding high-flow states have been shown to be sta- 
ble iQ (subject to the model's dynamics). In contrast in 
celltnar automata models as the VDR-model [^|l^] they 
usually show metastable behaviour. 

For application purposes open boundary conditions 
play a cruicial role. But also for purely theoretical pur- 
poses the investigation of the Krauss-model's behaviour 
under open boundary conditions is valuable. In general, 
driven interacting particle systems show boundary-induced 
phase transitions. For systems with unique flow-density 
relation, e.g. the ASEP as the prototype of such mod- 
els, there exists a quite general theory for the stationary 
state the system realized with open boundary conditions 
(cf. section ^) . These results are well established for sim- 
ple systems with unique flow-density relation [p7||2^ ,^. 
In contrast, not much is known about for systems with a 
non-unique fundamental diagram. Moreover, the models 
investigated so far have Wmax = 1 (i-e. Wmax can be reached 
within one timestep). Only recently also discrete models 
with higher velocities (e.g. models of Nagel-Schreckenberg 
type and non- unique flow-density relation p6|,p^ 

have been investigated. Here for the first time a continu- 
ous model with t^max 7^ 1 and fundamental diagram with 
bistability, i.e. stable high-flow states, has been studied 
with respect to boundary-induced phase transitions. 



In order to guarantee that the system remains free of 
crashes, one has to find a strategy to inject cars into the 
system. In this article two methods are discussed. The 
first one is orientated on the inflow rate a, i.e. with prob- 
ability a a car is always inserted if there is enough space 
at the left boundary. In this case one has to find an ini- 
tial velocity which guarantees a crash-free motion of each 
inserted car at any timestep (cf. section The second 
method fixes a minimum free space ginit at the beginning 
of the system ((7init > ^cai) and cars are injected with the 
constant velocity Umax (cf. section 

For both rules the phase diagram has been derived 
from computer simulations and, as expected, boundary- 
induced phase transitions were found. From the ASEP 
three different phases are known, distinguished by the 
functional dependence of the current through the system 
on a and f3. These are the low-density, the high-density 
and the maximum-curent phase. All these phases are ob- 
served in the Krauss-model with open boundary condi- 
tions. 

In contrast to the findings in the ASEP 0,|l|,^,|2^, 
^ the maximum current phase was only observed for an 
open right boundary (/3 = 1). This reflects the high sen- 
sitivity of the Krauss-model to an external disturbance 
at the system's exit which results from the long interac- 
tion horizon of the model. However, in contrast to |Q the 
maximum current phase in the Krauss-model exists for ar- 
bitrary system lengths due to the stability of the high-fiow 
states. 

Moreover, the extremal principle (0) [^,^ has been 
checked for both rules. For systems with Umax = 1 it was 
found that the selection of the system's state does not 
primarily depend on the parameters (a, /?), but on the re- 
sulting densities p^ and p+ at the left and right boundary, 
respectively. Furthermore, the selected state is completely 
determined by the fiow-density relation of the correspond- 
ing periodic system. In that sense the extremal princi- 
ple states implies the independence of the system state 
from the specific injection/extraction rule at the bound- 
aries. For the second rule we found absolut agreement for 
the Krauss-model (recall that here several timesteps are 
needed to accelerate to Umax)- 

However, for the first rule deviations were found from 
the principle. An unusual reentrance transition to a second 
low-density phase occurs for high insertion and extraction 
rates. Its existence could be ascribed to the behaviour of 
the rule which shows a cutoff in the maximum flow for 
a > 0.6. This flow is lower than the one that is allowed 
by the right boundary conditions. In other words, for a 
model with Umax 7^ 1 the selection of the system's state is 
not independent of the specific boundary rules. We there- 
fore argue that the extremal principle as given in [ p8| , [29| ] 
can only be applied in this strict form to models with 
t'max — 1- Here the provided maximum flow by the rule 
plays no role, since a stopped car can always accelerate to 
fmax in one timestep. Therefore any flow limitation due 
to the insertion rule does not play a role. This is different 
in the case of our model which has Umax > 1 and finite ac- 
celeration capability. However, if one takes an additional 
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rule into account, i.e. that the flow in the system can not 
exceed the flow aflowed by the boundary rules, one can 
also predict the system state with open boundary condi- 
tions for our first rule (cf. section ||) from the knowledge 
about the periodic system. 

We will close with a brief outlook. Due to our find- 
ings more work on boundary-induced phase transitions 
for models with Wmax > 1 seems to be necessary. The fo- 
cus should be on finding a compact formulation of the ex- 
tremal principle that is independent of Wmax- Since in the 
Krauss-model crash-free motion is not included as a rule 
per se (as, through hard-core exclusion, in the cellular au- 
tomata approaches) it would be helpful to develop a model 
formulation that does not show such safety problems un- 
der open boundary conditions. This would be important 
for applications in traffic flow simulation. To take into ac- 
count next-nearest-neighbor interactions as suggested by 
the findings of section y seems to be promising ||32| . 
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